MIXED ELLIPTIC AND HYPERBOLIC SYSTEMS FOR THE EINSTEIN EQUATIONS 

Yvonne Choquet-Bruhat and James W. York 

INTRODUCTION. 

We examine the Cauchy problem^ for General Relativity as the time history of the two 
fundamental forms of the geometry of a space like hypersurface, its metric g and its extrin- 
sic curvature K. By using a 3+1 decomposition of the Ricci tensor, we split^'^ the Einstein 
equations into "constraints", equations containing only g, K and their space derivatives, 
and evolution equations giving the time derivatives of g and K in terms of space derivatives 
of these quantities and also of the lapse and shift, i.e., the choice of the time lines. The 
constraints can be posed and solved as an elliptic system by known methods. However the 
equations of evolution in time for g and K are not, despite their form^ as canonical hamil- 
tonian equations, a mathematically well posed system; and they do not manifest directly 
the propagation of physical gravitational effects along the light cone. These equations, of 
course, contain gauge effects and cannot, therefore, yield directly a physical wave equation, 
that is, a hyperbolic system with suitable characteristics. (For a review of the relevant 
geometry, see, for example^.) 

In this paper we give two different methods for obtaining a nonlinear wave equation for 
the evolution of the extrinsic curvature K. The evolution of g is just its dragging by K 
along the axis orthogonal to the time slices. 

The first method corresponds to a choice of time slicing by fixation of the mean extrinsic 
curvature of the space slices. It leads to a mixed elliptic and hyperbolic system for which 
we prove local in time, global in space, existence theorems in the cases of compact or 
asymptotically euclidean space slices. 

The second method, which relies on a harmonic time slicing condition, leads to equations 
of motion equivalent to a first order symmetric hyperbolic system with only physically 
appropiate characteristics. We construct this system explicitly. Among the propagated 
quantities is, in effect, the Riemann curvature. The space coordinates and the shift are 
arbitrary; and, in this sense, the system is gauge invariant. 

Our gauge invariant non linear hyperbolic system, being exact and always on the physical 
light cone, is well suited for use in a number of problems which now confront gravity 
theorists. These include large scale computations of astrophysically significant processes 
(such as black hole collisions) that require efficient stable numerical integration, extraction 
of gravitational radiation with arbitrarily high accuracy from Cauchy data, gauge invariant 
perturbation and approximation methods, and posing boundary conditions compatibly 
with the causal structure of space time. 
n+1 DECOMPOSITION, ARBITRARY SHIFT. 

The n+1 decomposition of the curvature of a lorentzian manifold and the resulting decom- 
position of Einstein equations into constraints and evolution equations have been given 
in the case of zero shift by Lichnerowicz (1939) and in the case of an arbitrary shift 
by Choquet (Foures)- Bruhat (1956). Arnowitt, Deser and Misner (1962) have obtained 
this decomposition through a hamiltonian formalism and have deduced from it interesting 
physical consequences, for example, the definition of the gravitational mass. 
We recall this straightforward decomposition, using a coframe with time axis orthogonal 
to the space slices and introducing the notation of spatial Lie derivative like, for example. 
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Fischer and Marsden^. We consider on a manifold V = MxH a pseudo riemannian metric 
of forentzian signature which reads 

with (t e H, x*, i = l,.-,n are local coordinates on M) 

0^ = dt , 9' = dx' + P'dt 

With this choice of coframe, the dual frame has a time axis orthogonal to the space slices 
Mf = M X {t} while the space axis are tangent to them. The Pfaff derivatives da with 
respect to 6* " are 

do = -§1- P^di , di = ^ 
The structure constants of our coframe are such that 



Hence 



We denote by f*^ the Christoffel symbols of the (time dependent) space metric g = 
gijdx^dx^ , and by lo^^ the coefficients of the riemannian connection of g in the coframe 
We have 

_ /\r-i/)„ AT ,.,0 



OQ — IV lyo^v , ^Oi — ^iO 

a;0. = ^N-\^og^J - gh(rdj)^'') = \N-\f^g,^ - Cpg,^) 
where we define the symmetrization by f{ij)= fij + fji (no factor of 1/2). 
We define for any t dependent space tensor T another such tensor doT of the same type 
by setting 

^0 = it - ^13 
where jCj3 is the Lie derivative on Mt with respect to 
The extrinsic curvature of the space slices is given by 

K,j = -Nu% =- \N-^d^gij 

and we find 

ujl = -NKf , ujl^ = -NKf + di(3k 

Curvature tensor 

The components of the curvature tensor of g in the frame are found to be, using the 
general formulas given in Choquet-Bruhat and DeWitt-Morette I, p. 306, (with an opposite 
sign) 

RL = RAi - KiKu + KjK,, 
= -{N-^d^K.j + K.^Kf + N-^Vjd,N ) 

Ricci tensor. 

From these formulas result the following ones for the Ricci curvature, where we have 
denoted by H the mean extrinsic curvature of the space slices, i.e., we have set 
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H ^ Kl 

R.. = -N-^doKij + HKij - 2Ki^Kf - N'^VjOiN + Rij 
R,^^N{-\/HK^ + djH) 

i?oo = N{V%N - NKijK'^ + doH) 
SECOND ORDER EQUATION FOR K. 

In the formula above R^ is, like the right hand side giving its decomposition, a t dependent 
space tensor, the projection on space of the Ricci tensor of the space time metric. We 
compute its do derivative. First we compute d^Rij. The infinitesimal variation of the 
Ricci curvature corresponding to an infinitesimal 5g variation of the metric is 

5R,, = \{V^V(,5g^^h - ^''Sg^j - V,d,{g''^6ghk)} 
This expression applies to {d/dt)Rij with dgij = {d/dt)gij and to CpRij with 5gij = Cpgij. 
Therefore, using the relation between dogij and Kij , we obtain 

= -V^i{NV\Kj^h) - V^i{K^^dhN - 2NR\^^KJ^ - NR^^.Kf^ 

+ VhV^iNKij) + Vjdi{NH). 
We now use the expressions for Rq^ and Rij to obtain the identity 
^ij = doRij — V(i-Rj.)Q = 

-do{N-^doK,j) + do{HK,j - 2K,^Kf) - do{N-^W jd,N) - NVidjH 

-V(,(K,);,c>^iV) - 2NR%^KJ^ - NR^^.K^ + V/,V'^(iVKy ) + HVjdiN 
This identity shows that for a solution of the Einstein equations 



(E) RotP ^ Pai3 

the extrinsic curvature K satisfies a second order differential system which is quasidiagonal 
with principal part the wave operator, except for the terms VidjH. The other unknowns g 
and N appear at second order except for the term doVjOiN . Wc will use these facts in the 
following paragraphs to obtain geometrical well posed systems for the Einstein equations. 

I. MIXED HYPERBOLIC AND ELLIPTIC SYSTEM FOR g, K, N WHEN H IS GIVEN. 

A procedure to reduce the second order equation for K obtained above to a quasi diagonal 
system with principal part the wave operator is to replace in the term VjdiH the mean 
curvature H by an priori given function h: this is a gauge condition on the space slices. It 
was used by Christodoulou and Klainerman^, with h = 0, in the asymptotically euclidean 
case. With this replacement the equations Qij = Qij become, when N is known, a quasi 
diagonal second order system for K with principal part the wave operator, namely: 
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with 

uK.j = -do{N-^doK,,) + V'^VhiNK,,) 
where Pij depends only on K and its first derivatives, on N and 9oN together with their 
space derivatives of order < 2, and is given by: 

P,j = doi-HKij + 2g^^KihKj^)) + do{N-^VjdiN) 

+ V^iiKj^hd'^N) + 2NR\^KJ^ + NR^^^Kf^ - HVjdiN + NVjdih 

while Qij, zero in vacuum, is: 

Qij = dopij - V(j(pj)o) 
When N and the sources p are known, the above equation together with 

ig') doQij = -2NKij 

constitute a third order quasi diagonal system for g, hyperbolic if N'^ > and g is properly 
riemannian. 

On the other hand the equation Rq = Pq together with H = h imply the equation 

- {K.jK^^ - pg)iV = -doh 

This equation is an elliptic equation for N when g , K and p are known. 
Note that for energy sources satisfying the energy condition we have — Pq > as well as 
K.K = KijK'^^ > 0, an important property for the solution of the elliptic equation. The 
mixed hyperbolic elliptic system that we have constructed will determine the unknowns N 
and ^ in a neighbourhood of M in MxH when the shift /? is chosen. 

LOCAL EXISTENCE THEOREM(case of given H). 
We consider the vacuum case. The system for g, K 

(1) dogij = -2NKij 

(2) nKij = 

is equivalent to a third order hyperbolic system. It requires for its solution only local 
Sobolev spaces (cf. Leray^), but the lapse N is now determined by an elliptic equation on 
each Mt which must be solved globally: 

(3) V'diN - K.KN = -doh 

We will consider the two cases of a compact and of an asymptotically euclidean M. In both 
cases we will use an iteration scheme solving alternatively the elliptic equation and a linear 
hyperbolic system deduced in an obvious manner from (1), (2). 

We specify the results in the physical case of space dimensions n = 3 . There are problems 
in extending them to other dimensions in the asymptotically euclidean case. 
Case of a compact M. 

In the case of a compact M a good choice in order to have a positive lapse is to take 
doh = f > 0. Standard elliptic theory gives the following lemma, where Sobolev spaces Hg 
are relative to some smooth fixed metric e on M, I is an interval of IR and the shift P is 
taken smooth on Mxl. 
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Lemma 1. a) If g{t) G H^, K{t) e H21 with K.K ^ 0, and f G H2 the lapse equation has a 
unique solution N(t) G -ff4. By the maximum principle this solution is such that N > on 
M if f > on M and f ^ 0. 

b) If moreover ge n C^-^{I, Hk), while f , X G n C^-''{I, Hu) 
l<k<Z 0<k<2 

then 

ng n c2-fc(/,iy2+fe). 
< < 2 

Remark. We have K.K > hence K.K ^ if if ^ 0; we will prove later that the 

equality H = h holds for a solution of the coupled hyperbolic and elliptic system if it holds 
initially, but we cannot use this property in the iteration scheme. The condition K.K ^ 
must be deduced from such an hypothesis on the initial data and the continuity properties 
possessed by solutions of the hyperbolic system. 

We now consider the linear system for ^(n+1), ^(n+l) 

(In) doQijin + 1) = -2N(n)Kij{n), 

(2n) □„%(n+l) = Pi,-(n) 

obtained by replacing in □ the metric ^ by a metric ^(n) and in also K by K{n) and 
N by N(n). 

Lemma 2. Hvpothesis. Given: 

1) a) N(n) > No > , N(n) G n C^-''{I, H2+k), 

0<k<2 

b) hG n c^-''ii,Hk),i=[o,T] , 

2<k<3 

c) g{n)e n C^-''{I,Hk) and K{n) e n C'^-''{I, Hk). 

l<k<S 0<k<2 

2) Cauchy data such that 

gij{n + 1)(0, .) = 7y G H3, Kij{n + 1)(0, .) = % G H2, 

where the metric 7 is uniformly equivalent to the given metric e and /c./c ^ on M. 
The data doKij{n + 1)(0, .) = kij, determined by the equation -Rij(0, .) = 0, belong to Hi. 

Conclusion. The Cauchy problem for the system (In), (2n) has a unique solution 

^(n+l)G n C^-^{I,Hk), K(n+1)G n C^-\l,Hk). 
l</c<3 0</c<2 

Moreover there exists an interval Ii — [0, Ti], < Ti < T, and a number e > depending 
only on the norms of the given quantities in their respective spaces such that ^(n+1) is 
uniformly equivalent to 7 on Mx Ii and 

II K{n + l).K{n + 1) ||co(Mx/i) > e • 

Proof. The system is equivalent to a third order hyperbolic system, and the results can be 
proved by standard methods. 

Remark . When ^(n), i^(n), N(n) satisfy the indicated hypothesis the equation (2^) is a 
wave equation for K(n+1) whose coefficients satisfy the required hypothesis for the exis- 
tence of a solution in fl C^~'^(/, H^). The integration of equation (1„) determines 

< /c < 2 
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then ^(n+1), but the general theory of first order equations is not sufficient to assert its 
required regularity. One must use the fact that g{n + 1) satisfies a third order hyperbolic 
equation to prove that result. 

We remark also that the regularity hypothesis made on h is stronger that the regularity 
we prove on K. 

Theorem. Let (M,e) be a smooth compact riemannian manifold. Let be given on Mxl, 
with I = [0,T], a pure space smooth vector field /? and a function h such that 

he n C3-'=(/, Hk) , doh >0,doh^O . 

2<k<3 

There exists an interval J = [0/], £ < T such that the system (1), (2), (3) has one and 
only one solution on MxJ 

ge n C^-'^{J,Hk),Ke n C^-'^{J,Hk), N e n C^-'^{3,R2+k). 

l<k<3 l<k<2 0<k<2 

with N > and g uniformly equivalent to e, taking the initial data 

gij{0, .) = 7y e Hs, Kij{0, .) = kij e H2 

if 7 is a properly riemannian metric uniformly equivalent to e and k.k ^ 0. 

Proof . One determines N(0) on Mxl by the equation on each 

A ^N{0) - k.kN(0) = -doh(0) 

and then proceeds to the iterations, with N(n) determined from g{n), K{n). 
One uses the elliptic estimates for N(n) and the energy estimates of first and second order 
for K{n) and g{n) to prove the boundedness of the iterates in the indicated norms, as 
well as to show that there exists h = [0,£i] C I such that, on Mxli, g{n) is uniformly 
equivalent to e and || K{n).K{n)\\ co > e > 0. 

The elliptic estimates and the first order energy estimates applied to the difi^erence of 
two successive iterates show that there exists J = [0,£] such that the sequence converges 
strongly but in a weaker norm 

(C0( J, H2) n Ci( J, Hi)) X (C0( J, Hi) n C'{J, Ho)) X (C0( J, H3) n C'{J, H2)). 

Since the sequence is bounded in L°°{J,H^) x L°°{J,H2) x L°°{J,H^), the dual of a 
Banach space, there exists a subsequence which converges weakly in that space, and the 
limit found from the strong convergence belongs therefore to that space. An analogous 
argument applied to the time derivatives shows the existence of a subsequence whose 
first [respectively second] time derivatives converge weakly in L°°(J x H2) x L°°{J x Hi) 
[respectively in L°°(J, Hi) x L°°{ J x -f^o)], and the limit found belongs therefore to these 
spaces. It can be proved using properties of strong and weak convergence that the limit 
satisfies the non linear system (1), (2), (3). It can also be proved that this solution is in 
fact in the space given in the theorem (cf. a proof for second order systems by Hughes, 
Kato and Marsden^'^ and for symmetric hyperbolic first order systems by Majda^^). 

Case of (M,e) euclidean at infinity. 

We now consider the case of a manifold M which is the union of a compact set and a 
finite number of disjoint sets ("ends") diffeomorphic to the exterior of a ball in IR"^; M is 
endowed with a smooth properly riemannian metric e which reduces on each end to the 
euclidean metric. 
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We will look for an asymptotically minkowskian solution of (1), (2), (3) on Mxl, i.e. a 
solution belonging to some weighted Sobolev space Hg^s on each space slice M^. We recall 
that Hg^s is the completion of Cg" in the norm 

where o" ^ = l+d^, d the distance in the metric e to some fixed point in M. 

The theory of elliptic equations on an asymptotically euclidean manifold in Hg^s spaces 

(Choquet-Bruhat and Christodolou^^) gives the following lemma. 

Lemma 1. a) If g{t) — e e H^^-i, K{t) G -^2,0 and f = d^h e ^2,1 the lapse equation has 
a unique solution with 1 - N(t) e H4 _i. By the maximum principle this solution is such 
that N > on M if ^o/i > on M. 
b) If moreover: 

g-ee n C^-'^{I,Hk,-i), Ke n C^-'^{I,Hk,o), 
l<k<3 0<k<2 

fe n C^-''{l,Hk,i) 

0<k <2 
then 

1-Ne n C-^il.H^+k,-!). 

0<k<2 

Proof, a) The multiplication properties of weighted Sobolev spaces show that with the 
hypothesis made {K.K)(t) G -ff2,i- Since - | < -1 < - | the equation 

V'd^u - K.Ku = -doh + K.K , N = 1 + u 

has one and only one solution u G H4 _i; one deduces the property N > from the 
maximum priciple and the fact that N tends to 1 at infinity . 

b) The first time derivative of u can be proven to be in i^s^j, but again with — | < 5 < 
- |, hence in H^ _i (not in H^^), and analogously for the second time derivative with 3 
replaced by 2. 

Lemma 2. Hypothesis. We give on Mxl, I = [0,T]: 

a) /?G n C>'-^{l,Hk,-i), he n C^-^l^H^,^) 

l<k<Z 2<k<3 

b) N(n) > No > , N(n) G H C^->^{I, H2+k,-i) 

0<k<2 

c) g{n)-ee D C^-^{I,Hk,-i),K{n) e D C^-\l,Hk,o), 

l<k<3 0<k<2 
with moreover dtKn^ fl C^~'^{I, Hk,i) . 

0<k<l 

d) The Cauchy data (independent of n) such that 

gij{n + 1){0, .) = -y^J, 'y - e E Hs-i, Kij{n + 1){0, .) = kij e H2,o, 
and the metric 7 is uniformly equivalent to the given metric e. 

The data doKij{n+ 1)(0, .) = kij, determined by the equation i?y (0, .) = 0, belong to 
Conclusion. The Cauchy problem for the system (In), (2n) has a unique solution such that 
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g{n+l)-ee D C^->'{I,Hk,-i),K{n+ 1) e D C^-f^il, Hk,o) 
l<k<3 0<k<2 

dtK{n + l)e n C^-\l,Hk,i). 
< /c < 1 

There exists an interval Ii = [0, Ti], < Ti < T, and a number e > depending only on 
the norms of the given quantities in their respective spaces such that ^(n+1) is uniformly 
equivalent to 7 on Mx/i. 

Proof. The first weighted energy estimate needed to prove the lemma is obtained for 
Di^(n+1) (denoted K to simplify the writing) and dfK by multiplying its wave equation 
by cr^Voi^*"', after checking that under the hypothesis made on the n*'* iterates we have 
aPij{n) G 0*^(7, L^). The estimates for dfDK and D^iiT are obtained after derivation 
of the equation in the space directions. We obtain their boundedness in i/o,2 because 
a'^DPij e C°(/, L^).The estimate for d^^K is obtained only in Ho,i because only adtPij G 
C°(/, L^), due to the term in B'^d^N. 

The estimate for K can be obtained by integration in t from the estimate of dtK: it implies 
K-keC^{I, Ho,u hence K G C^{I, ifo,o) if it is so of k. 

The property dtg{n+l) G CO(I,H2,o), d lg{n+l) e C%I, d L^(n+1) G C^{I, ifo,i) 

results from the equation (In)- By integration with respect to t we obtain g{n + 1) — 7 
G C°(iy2,o)) hence g{n + 1) — e G C°(/, 1^2,-1) by the hypothesis on 7. To prove that it 
belongs to C°(/, -f^3,-i) one must use the third order equation for g{n + 1) deduced from 

(ln),(2n). 

The convergence of the iterates in a small enough time interval is proved analogously as 
in the case of a compact M and gives the following theorem. 

Theorem. The system (1), (2), (3) with Cauchy data 7 and k on the manifold (M,e), 
euclidean at infinity with given /3 and h on Mx [0,T] satisfying the hypothesis spelled out 
in lemmas 1 and 2, has one and only one solution (^,K,N) on MxJ, J = [0/] a sufficiently 
small subinterval of I, if the Cauchy data are such that 7 - e G H3 _i and is uniformly 
equivalent to e, k G H2,o, k G Hi^. The solution belongs to the functional spaces indicated 
for the iterates in the lemmas and is such that N > and g is uniformly equivalent to 7. 

Remark. The estimates show in fact that — 7 is such that for each t we have 

- 7 G i?3,o, 

i.e., it has a stronger asymptotic fall off than 7 - e or — e; this property is related to the 
A.D.M.^ theorem of mass conservation. 

SOLUTION OF THE FULL EINSTEIN EQUATIONS. 
We suppose that we have solved the reduced equations, that is 

nij + NVidj{H-h) = Qij 

and 

Roo - Ndo{H -h)= Poo 

We will show that the solution obtained satisfies the original Einstein equations if the 
source p is such that the associated stress energy tensor satisfies the conservation laws 
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VaT"^ = 

where (SttG = c =1) 

and the initial data 7 and k satisfy the constraints: 

50 = TO, Sf = T^, for t = 

with moreover 

trk = H| t=o = h| t=o 
while the initial values of OqK are chosen so that 

Rij |t=o= k=0) 
i.e., the Einstein equations are satisfied initially. 

The proof follows with Aaf3 = -Ra/3 — Pa/3) H — h = X. We will use the Bianchi identities 
Va<S"^ = where 

For any covariant symmetric space time 2-tensor A we find by using the values of the a;'s 
that we have the identities 

VoAij = doAj + N^Aj^h -N-^d^,NAj)o 
V^Ajo = ViAjo + N-^KijAoo - N-^Aojd,N + NAhjK'l 



A straightforward computation using these identities shows that the reduced equations can 
also be written: 

VoAij - V(jAj)o + NVidjX + KijdoX = 

and 

Aoo - NdoX = 

The Bianchi identities together with the conservation laws give the conservation equations: 

VaS«^ = 0, with S"'^ = - 

In our frame we have the identities 

-^ij = + yi^{N'^A°° - g^^Ahk) 
The conservation equation with = implies therefore 

i(Vo^°° + N-'^g'^VoAij) + Vi^^o = 0. 
Hence, when the reduced equations are satisfied (recall that VqN = 0), 

-N-'^VodoX + g'^VidjX + N-^HdoX = 0. 
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This wave equation for X implies X = OonMxIifX and OqX vanish initiaUy. We have 
then also Aqq = 

We can now proceed in a manner analogous to ^. We differentiate the conservation equation 
with (3 = i and use the Ricci identity to obtain 

VoVo^o^' + V.VoA^'^ + Ro,\A^' - RoxA^^ - ^V^VoA\ + ^N^V^VoA'"^= 0. 

Using the reduced equations now with Aqq = X = 0, we find 

(VoVo - N'^V^V')A^^ - iV2(v,VM*o - \/^V^A'^) + iA^2yjy^^00 _ Q 

where ~ means modulo linear terms in A^^. Using again the Ricci identity and the 
conservation law with /? = 0, we find 

-□ A^'o = (TV-^VoVo - ViV')A^^ ~ 

where ~ denotes terms linear in and its derivatives of order < 1. We deduce then 
from the first reduced equation: 

□Vo^y ~ 

where ~ is modulo linear terms in A'^^ and their derivatives of order < 2. We have 
obtained for Aix a strictly hyperbolic system which is linear and homogeneous. The 
uniqueness theorem of Leray^ for solutions of such systems shows that we have Aix = 
on Mxl if the corresponding initial data are zero. 

The set of initial data in our problem are the values on Mq of X, OqX, A^q, dgAiQ , 
Aij, doAij, dggAij. These values can be proved to be all zero for a solution of the reduced 
equations if {H — h) \t= and the Einstein equations are satisfied on Mq. Our proof is 
complete. 

II. HARMONIC TIME, ALGEBRAIC GAUGE 
HYPERBOLIC SYSTEM FOR g, K, N. 

We shall eliminate at the same time the third derivatives of N and the second derivatives of 
H in the second order equation for X by a particular choice of N which we now explain (cf. 
this elimination in the case of zero shift in Choquet Bruhat and Ruggeri^). We compute 

doVjdiN={^^-Cp)VjdiN 

We find at once 

and we find an anlogous formula when the operator d/dt is replaced by >C/3 by using the 
property 

diCi3N = CpdiN = 

which leads to 

= -\g^\V^iCpgj)i - ViCpgij)dkN 

Finally 
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The third order terms in N and second order terms in H in the second order equation for 
K can therefore be written under the following form 

dj = -N-^Vjdi{doN + N^H) 

We satisfy the condition = by imposing on N the differential equation 

(N') doN + N'^H = 

The second order equation for K reads then as a wave type non linear system: 

(K') NnKij = NQij + Qij 

where we now set 

nKij = -N-^dodoKij + V^VhKij 

NQij = -KijdoH + 2g'''^KmiidoKj)h + 4.Ng^^g'^^KikKimKjh 
+ {2V^iKj)i)d^N - 2HN-'^diNdjN - 2d^iNdj)H 
-ZdhNW^'K.j - K,,W^WhN - N-^Kijd^NdhN + N'^ K^^d^^Nd^N 
+ (V(ia'^iV)X,)ft + 2NR%^K]:' + NR^^.Kf^ - 2HV ,d,N 

Sij = dopij - ^{iP^^) 

The equation (N') expresses that the time coordinate is harmonic, namely: 

V^d^x^ = N'^oj^Q - g^^uj% = 0. 

Using the expression for H we see that this equation reads 

dQlog{N/{detgf'^} = 

We find, generalizing the result obtained by Choquet-Bruhat and Ruggeri^ in the case of 
zero shift, that the general solution of this equation is 

N = a-\detg^/^) 

where a is an arbitrary tensor density such that 

doa = 0. 

This equation is a linear first order partial differential equation for a on space time, de- 
pending only on (3. A possible choice if P does not depend on t is to take a independent 
of t and such that £^30; = 0. 
The above choice of N is called algebraic gauge. 

If we replace N by the value obtained above in the second order equation for K, we obtain 
a quasi diagonal system with principal part the wave operator, with terms depending on 
g and its derivatives of order < 2. This system reduces to a third order hyperbolic system 
when we replacing K by —{2N)~^dog. The local existence theorem of Leray^ for the 
solution of hyperbolic systems gives immediately the local in time existence of solutions of 
this reduced system, in local in space Sobolev spaces. 

It can be proved that a solution of the reduced system on Mxl is a solution, in algebraic 
gauge, of the full Einstein equations if the initial data satisfy the Einstein equations on 
Mo (cf. the proof in the case of zero shift in^). 

FIRST ORDER SYSTEM (vacuum). 
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The preceding results can be extended without major change to dimensions greater than 4. 
We will show that in dimension 4 a solution of the vacuum Einstein equations, together with 
the harmonic time gauge condition, satisfies a first order symmetric system, hyperbolic if 
g is properly riemannian and > 0. Such a system could be useful to establish a priori 
estimates relevant to global problems. It is of great importance for numerical computations, 
because symmetric hyperbolic systems occur in many areas of mathematical physics, in 
particular in fiuid dynamics, and codes have been developed to study such systems. 

We have obtained for the unknowns g, K, N the equations 

(1) dogij = -2NK,j 

(2) doN = -N^H 

(3) □ Kij = Qij 

To obtain a first order system we take as additional unknowns: 

doK,, = Mq,j, VhK^j = Mhij, diN = N^, BodiN = Noi, VjdiN = Nji. 
We take as equation (3') 

(3') doKij = Moij 

The equation (3) gives 

(4) doMoij - N^V^Mhij = -N^Qij 

In three space dimensions the riemann tensor is a linear function of the Ricci tensor: 

Rlijm = dlj-^im ~l~ 9imRjl dij-^lm 9lmRij '^{,9lj9im 9ij9lm)R 

Using the equation i?y = we have 

We then can write N'^Qij as a polynomial in the unknowns, and g^^ . 

We prove the following lemma for instance by using the commutativity dodi = dido on 
components of tensors and the value of dof^j in terms of K: 

Lemma. For an arbitrary covariant vector u^ we have 

doVhUi = Whdoui + ui{Wh{NKl) + Vi{NK{) - V\NKih)} 

and an analogous formula for tensors with additional terms for each index. 

Using this lemma we see that Moy and M^ij must satisfy the equations 
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On the other hand by definition 

(6) doNi = Noi 
while Nhi and Noi must satisfy 

(7) doNhi - VhNoi = NNiiMj + - M\^) + Ni{kINh + K^Ni - N^K^h). 

Now we deduce the value of dodoNi = doNoi from (2) and from the Einstein equation 

i?o = -{N-^V'^VhN - K,jK'^ + N-^doH} = 

Indeed these two equations imply that N satisfies the following inhomogeneous wave equa- 
tion 

Hence by differentiation, use of the Ricci formula and the definitions of Nhi and Noi, we 
find 

(8) 

doNoi - N^V'^Nm = -R'^NnN^ + 2NN,N\ - 2N^KmM^' 
-SN^NiKhiK''^ + AN^HM^X + QN^H^Ni. 

We use again the equation Rij = to replace Rij by its value in terms of the unknowns. 

We have obtained a first order system (equations 1, 2, 3', 4, 5, 6, 7, 8) in all the unknowns. 
The right hand sides are polynomial in the unknowns, g'^^ and N~^. They do not depend 
on their derivatives. The left hand sides are linear operators on the all the unknowns. 
Their coefficients depend on these unknowns, and not on their derivatives except for the 
derivatives of g. In order to obtain a covariant quasilinear first order system for all the 
unknowns we can, for instance, introduce on M an a priori given metric e which may 
depend on t but is such that 

doCij = 0. 

We denote by D the covariant derivative in the metric e. We introduce the additional 
unknown G^ij = D^gij and deduce from (1) the equation: 

(9) doGmj = -2{NhKij + N{Mhij + S^^.K^)^)} 

where the tensor S, the difference of the connections of g and e, is given by 
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We have now obtained a quasi linear first order system for all the unknowns. Its charac- 
teristic matrix, obtained by replacing in the principal matrix the operator 9 by a covariant 
vector ^, is constituted of blocks around the diagonal, some reduced to one element ^o, and 
some 4x4 matrices with determinant — N'^^^^i. The characteristics are the light cone 
and the time axis. On the other hand the system can be symmetrized by multiplication 
with a matrix constituted of blocks around the diagonal equal to one element, 1, or the 
matrix [g'^^). In other words the system is a first order symmetrizable hvperbolic system, 
with domain of dependence determined by the light cone . Known local existence theorems 
apply to such a system. 

Remark. It is possible to obtain a system whose right hand sides are polynomial in the 
unknowns and g'^^ by making a slightly different choice of these unknowns, in particular 
by introducing = N-'^diN (cf.^3,i4^_ 
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